Next to leading order non Fermi liquid corrections to the neutrino emissivity and 

cooling of the neutron star 

Souvik Priyam AdhyaQ Pradip K. RoyQ and Abhee K. Dutt-MazumdeiH 

High Energy Nuclear and Particle Physics Division, 
Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata-700 064, INDIA 

In this work we derive the expressions of the neutrino mean free path(MFP) and emissivity with 
non Fermi liquid corrections up to next to leading order(NLO) in degenerate quark matter. The 
calculation has been performed both for the absorption and scattering processes. Subsequently 
the role of these NLO corrections on the cooling of the neutron star has been demonstrated. The 
cooling curve shows moderate enhancement compared to the leading order(LO) non-Fermi liquid 
result. Although the overall correction to the MFP and emissivity are larger compared to the free 
Fermi gas, the cooling behavior does not alter significantly. 
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I. INTRODUCTION 

It is known for quite sometime now that the degenerate Fermi gas at fow or zero temperature gives rise to phe- 
nomenon which is different from the normal Fermi hquid behaviour once the magnetic interactions are included. This 
interesting feature is characterized by the appearance of anomalous logarithmic term in the expressions for various 
physical quantities like specific heat, entropy etc[l|. It has also been revealed recently that at strictly zero tempera- 
ture, there is a logarithmic singularity in the inverse group velocity, which leads to the breakdown of the usual Fermi 
liquid picture in presence of the magnetic interactions. Historically, such a deviation from the normal Fermi liquid 
behavior was for the first time, exposed in where, the specific heat of a degenerate gas due to the current-current 
interactions was calculated and the result contained the T\iiT~^ term which emanates from the unscreened magnetic 
interactions. It is to be mentioned here that non-Fermi liquid behaviour of highly dense color superconducting QCD 
plasma has been studied in great detail [I, Q . 

For non-relativistic systems, the magnetic interaction is suppressed in powers of (vjcf' and therefore might not 
be of much quantitative importance. However, for dense plasma where the constituents like quarks or electrons are 
moving with a velocity close to the velocity of light, the magnetic interactions cannot be neglected. In fact, it has 
been revealed recently that in many context, the transverse interactions, due to its infrared sensitivity, may become 
more important than its longitudinal counterpart in this kinematic regime. For example, while calculating the fermion 
damping rate and energy loss, it has been shown in ref.0, Q that the leading order (LO) contributions come from the 
magnetic interaction while the longitudinal interactions contribute only at the sub leading order. In fact, it has been 
seen that the first two leading order contributions in the expressions for the fermion dampingrate in ultradegenerate 
plasma come from the transverse sector alone. Similar behaviour has also been reported in where the authors 
have studied the non- Fermi liquid behaviour (NFL) of the drag and diffusion coefficient in degenerate plasma. A more 
elaborate discussion on the NFL aspects of the cold and dense QED and QCD plasma has been presented in 

This recently discovered phenomenon of non- Fermi liquid behaviour, which relates itself to the modified quark 
dispersion relation for excitations close to the Fermi surface, also finds important application in astrophysics. For 
example, it has been shown that the NFL corrections to the quark self-energy enhance the neutrino emissivity of 
ungapped quark matter which may exist in the core of neutron stars0, [T3|. Like emissivity, in dense quark matter, 
the neutrino mean free path(MFP) also receives significant NFL corrections as has been demonstrated in It might 
be mentioned here that in all these calculations the evaluation of the quark self-energy was restricted to the leading 
logarithmic order. In [T2j . on the other hand, the authors determine the quark dispersion relations in ultradegenerate 
relativistic plasmas beyond LO, which, at zero temperature, is characterized by the appearance of the fractional higher 
powers in the energy variable. In another work, the specific heat of normal degenerate quark matter has also been 
calculated where also in the higher order terms the fractional powers show up [l3| . 

In view of these contemporary investigations, we here plan to evaluate the neutrino MFP and corresponding 
emissivity in normal degenerate quark matter beyond leading logarithmic order and compare with the LO results. 
Here leading order (LO) refers to the anomalous logarithmic term Tlog(f/T) that occurs as the first term in the 
non-Fermi liquid contribution to the fermion self energy. Quantities such as mean free path and emissivity calculated 
with this term is called the LO corrections. Next to leading order (NLO) terms include all other terms beyond the 
LO that contain the fractional powers of T and up to the (T'^)log(I/r) that occur in the expression of the fermion 
self energy. Similarly, quantities calculated with this correction are labeled as NLO corrections [l2| . Equipped with 
these results and knowing the specific heat of dense quark matter upto the order concerned, we investigate the cooling 
behaviour of the neutron star with dense quark core. 

The plan of the paper is as follows. In section II, we describe the formalism; where we start with the quark 
dispersion relations and the modifications due to NFL effects followed by evaluation of the MFP for the degenerate 
and non-degenerate neutrinos. These are followed by the calculation of emissivity of the neutrinos. Section HI is 
devoted to the study of the cooling process via neutrino emission. Finally, the results are summarized in section IV 
followed by conclusion in section V. 



II. FORMALISM 

A. Quark dispersion relation 

To calculate the quark dispersion relation in degenerate plasma one needs to evaluate the quark self-energy. For 
this, we consider Fig.(IT]), where the solid line represents the fermion prop agator and the blob implies that the gluon 
propagator used here is hard dense loop (HDL) corrected propagator [Ijl . Mathematically, the quark self energy can 
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FIG. 1. Fermion self-energy with resummed gluon propagator. 



be written as d i, [H [H [H , 

d3 



(1) 



where, P'' = (p^jp), Pq = iuJn + /^j 9o = "i^s, w„ = (2n + l)7rT and u/g = SttsT are the Matsubara frequencies 
for fermion and boson respectively with integers n and s. After performing the sum over Matsubara frequency in 
Eq.(IT]), iuj„ + /i is analytically continued to the Minkowski space so that = (£^,p)[i2[. S(P) can be written as a 
combination of quasiparticle and antiquasiparticle self energies as[T2|. 



I](P) = 7oA+S+(P) - 7oA-S_(P) 



(2) 



where 

5±i = -[/t(|p| + s±)] 

(3) 

with the energy projection operators are given by, 

A± = l(l±7o7y) (4) 

Interactions within the medium severely modify the on-shell self-energy of the quarks which is manifested in the 
slope of the dispersion relation for the relativistic degenerate plasma. For quasi-particles with momenta close to the 
Fermi momentum pf(i) where i denotes the quark flavor, the one- loop self-energy is dominated by the soft gluon 
exchanges @ . For the calculation of the MFP and emissivity one needs to know the modified dispersion relation 
which is determined as 

w± = ± (Pp(„^ ) + RcE± {uj±,p{uj±))) (5) 

where w is the quasiparticle/antiquasiparticle energy which is a solution of the dispersion relation and £'p(w) = 

yjp'^luj) + is the kinetic energy. As we are considering only quasiparticles, we will consider only and denote 

a;-|_ by uj. The above expression will be used to obtain dp/du needed for the phase space evaluation of the mean free 
path. The authors of have already calculated the fermion self energy with terms beyond LO. We here quote the 
low temperature expansion of the on-shell fermion self energy for |p| = E (ultrarelativistic case) 0, [l^ and notice 
that no explicit dependence on the spatial momentum p occurs (isj : 

= ^ - ^^-^ { 12^ (^) + 1] + ^ + (-) (^g-(-) - ^^^) 



64a/2 Vm/ 1897rii/3 
6144 - 2567r2 + SGtt^ - Ott*' / e 



8647r6 



11/3 ~ 
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where e = {uj — n) ^ T where NFL effects dominate. The scale of the last logarithm was determined by resumming 
infrared enhanced contributions. The first term in the above expression gives the hard part contribution to the self 
energy where = g^Ci? /x^/(47r^) and m is given by = Nfg'^^'^/^Air^) and is related to the Debye mass by 
= m|)/2. It is interesting to note here that at higher order, fractional powers in e appear. This can be attributed 
to the dynamical screening for the transverse exchange of gauge bosons. 



B. Mean free path (MFP) of the degenerate neutrinos 



The degenerate neutrinos refers to the case where the neutrino chemical potential (/i^) is much larger than the 
temperature. In the interior of a neutron star, there are two distinct phenomena for which the neutrino mean free 
path is calculated, one is absorption and the other involves scattering of neutrinos [l3 ■ To calculate the MFP for the 
absorption process we consider the simplest /3 decay reactions; i.e the absorption process and its inverse p7l - [l9| . 



(7) 



(8) 



The corresponding mean free paths are denoted by l^^an ^^'^ ^mean- 

The neutrino MFP is related to the total interaction rate due to neutrino emission averaged over the initial quark 
spins and summed over the final state phase space and spins. It is given bv[l7l|. 
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(9) 



where, g' is the spin and color degeneracy, taken to be 6 and p and Up are the energy, momentum and distri- 
bution function for the corresponding particle. |Mp is the squared invariant amplitude and is given by |Afp = 
64G^ cos^ Oc(Pd ■ Pu){Pu ■ Pe)- Here G ~ 1.435 x lO^^^erg — cm? is the weak coupling constant. Here, we work with 
the two flavor system as the interaction involving strange quark is Cabibbo suppressed. We now consider the case of 
degenerate neutrinos i.e. when ^i^, ^ T. So in this case both the direct Eq.© and inverse Eq.® processes can occur. 
Consequently, the /3 equilibrium condition becomes + Mi^ = + /^e- Now, to carry out the momentum integration, 
cPpd and d?Pu can be evaluated as. 



,3 o Pf{d) , dpd 3 
a pd = 2n — —pdp——duj;d p^, 
Pf{v) did 



p duj 



(10) 



where we define p : 

El, 



\Pd +Pu\ — \Pu+Pe\- dp{u!)/duj can be evaluated from the modified dispersion relation as follows 

ai?eE+(w) doj 



du} 



dE, 



'p(") 



dp[uj) dp[uj) 
dp{uj) 



dio 



= 1 



doj dp(uj) ' 
dReT,+ {LLi)\ Ep(^^) 
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p{io) 
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where 9i?el]+/9p ~ 0, since p does not appear explicitly in the expression for Neglecting the quark-quark 

interactions, the leading order result is obtained as. 
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The NLO result is evaluated as. 
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where the constants are. 
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and 
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To arrive at the Fermi-liquid result, one can use the free dispersion relation to arrive at, 



(15) 
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Since quarks and electrons are assumed to be massless, the chemical equilibrium condition gives Pf{u) +pf{e) = 
Pf{d) +pf{i'), which we use to derive Eg. (fT6|) . We have further assumed that fid ^ jJ^u = A*- Next we calculate the 
MFP for the quark-neutrino scattering process. 



+ Ve{Ve) qi+ Ve{Ve) 



(17) 



for each quark component of flavor i{— u or d). Including the NFL corrections through the phase space and assuming 
raqjpf^ <C 1, we obtain. 
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where Uq^ is the number density of quark of flavor i, given by, 

d^p 1 



(27r)3 e'3(^«-A',. 



(20) 



(21) 



where 6 is the quark degeneracy factor, niq^ is the mass of quark, gq = AG'^ml/Tr 0, [2l| and A{xi) is defined in 
[ll| where Xi = fi^/ fiq. if fi^ < fiq. and Xi = fiqj fii, if /ij/ > /Xg. . The contributions from the Fermi liquid(FL), LO 
and NLO are added to obtain the MFP for the corresponding process. Further, one can combine l^tmi with l':^lan to 
define total mean free path as [22j . 

111 
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In addition to the inclusion of the NLO terms, the Fermi liquid term and LO term agrees with that of |ll| . 

C. Mean free path of non-degenerate neutrinos 

We now derive MFP for nondegenerate neutrinos i.e. when fi^ <^ T beyond the Fermi-liquid contribution. For 
nondegenerate neutrinos the inverse process (jS]) is dropped. We are considering only depopulation of neutrinos as in 
non-degenerate (untrapped) case repopulation or the reverse reaction is assumed to be zero. Hence, we neglect the 
second term in the curly braces of Eq. ([5]) [13, [l^l ■ For free quarks, the matrix element vanishes [13, El, 113' since u, 
d quarks and electrons are collinear in momenta. The inclusion of strong interactions between quarks relaxes these 
kinematic restrictions resulting in a nonvanishing squared matrix amplitude. We can neglect the neutrino momentum 
in energy conserving relation due to the thermal production of the neutrinos [l7j . Following the procedure described 
in [H, I25I, [2^ the MFP for the Fermi liquid case, LO and NLO are obtained as, 
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where the constants are evaluated as, 



and 



bi = 



6. = 



25/3 



280 X 2^/3 



9V37r5/3' " 189\/37r7/3 

6144 - 2567r2 + SGTr* - Ott^ 
2167r4 ■ 



(25) 



(26) 



(27) 



Similarly, for the scattering of nondegencrate neutrinos in quark matter with appropriate phase space corrections we 
obtain. 
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where the constants are. 
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Here, we have assumed mq./pf. <^ 1. Thus, the total MFP for non-degenerate neutrinos is obtained by summing up 
the contributions from the absorption and scattering parts to get the expression of the MFP of the non-degenerate 
neutrinos up to the NLO terms. So, we get the expression of the MFP of the non-degenerate neutrinos up to the 
NLO terms. 



D. Emissivity of non-degenerate neutrinos 

The total emissivity of the non-degenerate neutrinos is obtained by multiplying the neutrino energy with the inverse 
of the MFP with appropriate factors and integrated over the neutrino momentum. The relation between neutrino 
emissivity and the neutrino mean free path is thus obtained as |27| . 



J (27r)3 •'1{~E,,[ 



Using the mean free path for the non-degenerate neutrinos we obtain, 

e — sq ~ elo + £nlo 

where, 
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(35) 



is the usual Fermi liquid contribution which agrees with the result presented in ref.[l7|. At the LO we have obtained, 
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FIG. 2. Mean free path of degenerate neutrinos. The left panel shows a comparison between the Fermi liquid result and NLO 
corrections for the non-Fermi liquid effects. The right panel shows the reduction of the MFP due to NLO corrections. 



which is in agreement with the result quoted in ref.Q. Now, following the procedure in 0, [Tl|, we obtain the NLO 
contribution to the neutrino emissivity as, 
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0.656g/i 
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where the constants are evaluated as. 
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The NFL correction only appear in the phase space integral of the MFP [Tl| and subsequently in the expression of the 
emissivity 0. It is actually related to the unscreened transverse interaction The factor of can be understood 
easily. Naively, we can see that one power of T is obtained from phase space integral of a degenerate fermion. Further, 
we obtain a from the phase space integral of the neutrino. One power of T from the energy conserving S function 
is cancelled by a power from the emitted neutrino energy. It is to be noted that for excitations near the Fermi surface, 
the angular integrals give no temperature dependence. Now putting the standard values [ll| for the parameters, 
the corrections can be compared with the results given in 



III. COOLING PROCESS VIA NEUTRINO EMISSION 



The temperature of the neutron star with a quark matter core shows a dependency with time. To analyse the 
cooling of the star [28l - [30j . the specific heat capacity of the quark matter core needs to be taken into consideration 
along with the emissivity via the cooling process 17. |31|. 



du _ dudT _ dT _ _ 

dt ~ dT dt """"^^^i ~ ^ 



(40) 



where u is the internal energy, t is time and we have assumed that there is no surface emission. The NFL effects 
on the specific heat capacity of the degenerate quark matter has been calculated recently in [ll] . The specific heat 




FIG. 3. Same as FigOfor non-degenerate neutrinos. 



capacity for a non-color-superconducting degenerate quark matter is given as [13|, 
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and A^/ is the number of quark flavors. The contribution from the free part is given as, 



(41) 



(42) 



(43) 



Using the above expression for specific heat and the emissivity expression upto NLO, we analyze the cooling behavior 
of the neutron star. 



IV. RESULTS 



An estimation of the MFP of neutrinos with the temperature has been presented in this section. For this purpose, 
we have assumed the quark chemical potential to be 500 MeV. This is in well agreement with the high density ~ Gpo 
(po being nuclear matter saturation density) at the core. The chemical potential may well have a time dependence. 
This would gain significance for terms beyond order since T/{giJL) « 1. We have also taken /i^ = 15MeU and 
Us = 0.1. In left panel of Fig.® we note there is a considerable decrease in MFP of degenerate neutrinos due to NLO 
corrections over the Fermi liquid result. In right panel of Fig.([2|) LO correction is compared with the NLO correction 
and it is seen that the MFP with NLO correction is reduced marginally as compared to LO correction. Similar 
features have been observed in case of the MFP of non-degenerate neutrinos as displayed in Fig. ([3]). This marginal 
difference in the MFP between NFL LO and NLO corrections also leads to marginal difference in the emissivity for 
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FIG. 4. The left panel shows the emissivity of the neutrinos with temperature in degenerate quark matter. The right panel 
shows the behavior of the specific heat of the degenerate quark matter with temperature (Tg in units of 10^ K). 
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FIG. 5. The cooling behavior of neutron star with core as neutron matter and degenerate quark matter. The dotted line 
represents the Fermi liquid result, the solid line represents the non-Fermi NLO correction. The dash-dotted line gives the 
cooling behavior of the neutron star core made up of purely neutron matter. 



the two cases. This is shown in Fig.(j3]). These small reductions are reflected in the marginally enhanced emissivity of 
the non-dcgencratc neutrinos which has been shown in left panel of Fig. (21). We find that there is a modest increase 
in the emissivity of the neutrinos. The right panel of Fig.Q gives a comparison of the NFL corrections to the specific 
heat already reported in The complicated cooling equation cannot be solved analytically and we have resorted 
to numerical calculation. We observe that the cooling of neutron star is marginally faster in case of NFL (NLO) as 
compared to the Fermi liquid result (shown in Fig.®). 
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V. SUMMARY AND DISCUSSIONS 

In this work, wc have calculated the MFP of degenerate and non-degenerate neutrinos both for the scattering and 
absorption processes. We then find the expression for neutrino emissivity for non-degenerate neutrinos with NLO 
corrections. It is seen that both MFP and emissivity contain terms at the higher order which involve fractional powers 
in (T//i). We have found that there is a decrease in the MFP due to NLO corrections. We reconfirm that the leading 
order correction to the quantities like MFP or emissivity are significant compared to the Fermi liquid results. The 
NLO corrections, which we derive here, have however been found to be numerically close to the LO results. We 
have also examined the cooling behavior of a neutron star by incorporating NLO correction to the specific heat and 
emissivity which affect the results considerably compared to the simple Fermi liquid case. 
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